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Abstract 



We investigate a peculiar supersymmetry of the pairs of reflectionless quantum systems de- 
scribed by n-gap soliton potentials of a general form that depends on n scaling and n translation 
(— I ■ parameters. We show that if all the discrete energy levels of the subsystems arc different, the 

^ f superalgebra, being insensitive to translation parameters, is generated by two supercharges of 

differential order 2n, two supercharges of order 2n + l, and two bosonic integrals of order 2n+l 
composed from Lax integrals of the partners. The exotic supersymmetry undergoes a reduction 
when r discrete energy levels of one subsystem coincide with any r discrete levels of the partner, 
the total order of the two independent intertwining generators reduces then to 4n — 2r + 1, and 
the nonlinear supcralgebraic structure acquires a dependence on r relative translations. For a 
complete pairwise coincidence of the scaling parameters which control the energies of the bound 
^ ' states and the transmission scattering amplitudes, the emerging isospectrality is detected by a 

. transmutation of one of the Lax integrals into a bosonic central charge. Within the isospectral 

class, we reveal a special case giving a new family of finite-gap first order Bogoliubov-de Gennes 
systems related to the AKNS intcgrable hierarchy. 



1 Introduction 



Solitons and related topologically nontrivial objects such as kinks, instantons, vortices, monopoles 
and domain walls play an important role in diverse areas of physics, engineering and biology [U 
^ ! Ells]. Darboux and Backlund transformations, with their origin in the theory of the linear Sturm- 

■ Liouville problem and classical differential geometry, proved to be very effective in their study 

[His]. Darboux transformations on the other hand, underlie the construction of supersymmetric 
quantum mechanics [SI 17]. Via the Bogomolny bound and the associated first order Bogomolny- 
Prasad-Sommerfield equations [319], supersymmetry, in turn, turns out to be closely related with 
the topological solitons [TO] . 

Solitons and their periodic analogs appear as solutions of classical nonlinear integrable field 
equations, and by means of Lax representation [11] are related with finite-gap quantum systems 
[121 113|. As the latter class of the systems is characterized by nontrivial, higher derivative integrals 
of motion, one could expect that supersymmetric extensions of them should possess some peculiar 
properties. This is indeed the case [HI HSl HH 113 UHl IH] , and exotic supersymmetric structures of 
finite-gap systems found recently some interesting physical applications [201 Ell 1221 l23l 123] . 

The most known example of non-periodic finite-gap systems is given by a hierarchy of reflec- 
tionless Poschl- Teller potentials. The Schrodinger Hamiltonian with one-, two-, and, in general, 
n-gap Poschl- Teller potentials controls, particularly, the stability of kinks in sine-Gordon, ip'^ and 
other exotic (l-hl)-dimensional field theory models [fl El 123 1211 123 123 129] • These reflectionless 
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systems also appear in Gross-Neveu model [30\ I31|. The indicated hierarchy represents, however, 
only a very restricted case of a general family of n-soliton potentials. The latter corresponds to 
2n-parametric solutions of the Korteweg-de Vries (KdV) equation [21 HI [32] . 

In the present paper we study the exotic supersymmetry that appears in the pairs of reflec- 
tionless systems described by n-soliton potentials of the most general form. This will allow us to 
reveal and investigate new features of the peculiar supersymmetric structure of finite-gap systems 
that experiences radical changes in dependence on relation between the two sets of the scaling 
and translation parameters of the partner potentials. Analyzing different faces of supersymme- 
try restructuring, we detect a special family of supersymmetric n-soliton partner potentials which 
provides us with new finite-gap solutions for the Gross-Neveu model by means of the first order 
Bogoliubov-de Gennes system. 

The paper is organized as follows. In the next Section, we review the general construction of 
soliton potentials with the help of Crum-Darboux transformations, summarize the basic properties 
of the corresponding refiectionless quantum systems, and formulate precisely the problems related 
to supersymmetry of finite-gap soliton systems to be studied here. Section 3 is devoted to the 
analysis of supersymmetry of non-isospectral pairs of refiectionless one-gap systems with different 
bound state energy levels. In Section 4 we investigate the changes this supersymmetric structure 
undergoes in the isospectral case. Section 5 generalizes the results of Section 3 for the case of n- 
gap pairs with completely broken isospectrality. To clarify the supersymmetry picture in extended 
n-gap systems with partially broken and exact isospectralities, we study in detail the case of n = 2 
in Section 6. In Section 6.1 we review the properties of the generic n = 2 refiectionless systems 
to identify the ingredients to be important for further analysis. Then, in Section 6.2, we discuss a 
generalization of Crum-Darboux transformations that is related to alternative factorizations of the 
basic Crum-Darboux generators of order n > 1. The results of Sections 6.1 and 6.2 are employed 
in Sections 6.3 and 6.4 for analysis of supersymmetry in extended n = 2 systems with partial 
isospectrality breaking. Finally, in Sections 6.5, 6.6 and 6.7 we investigate the most tricky case 
of supersymmetry in two-soliton extended systems with exact isospectrality. We do this first in 
Section 6.5 for a particular case of exact isospectrality with a common virtual n = 1 subsystem. In 
Section 6.6 we investigate a generic case of exact isospectrality, within which we detect yet another, 
very special, particular case. The latter is studied in Section 6.7, and provides us with a new, 
first order finite-gap system belonging to the AKNS hierarchy [33[ 113] . In Section 7 we discuss 
how the results on partially broken and exact isospectralities are generalized for extended systems 
with n > 2. We conclude the paper with discussion of the obtained results and their possible 
developments and applications in Section 8. 

2 Family of refiectionless n-gap systems 

A Crum-Darboux transformation of order n, n = 1,2,..., applied to a quantum free particle 
generates a system characterized by the Hamiltonian [1] 



Here Hq = —-^ is a free particle Hamiltonian, and = VF(V'i) ■ ■ ■ i V'n) is a Wronskian of its 
(formal) eigenfunctions iIji{x), . . ., ipnix), HqiPj = Ej^^ipj, 








i,j = 1, . . . 



n. 



(2.2) 
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If the set i^jix), j = 1, . . . ,n, is chosen in a way such that ^ 0, the potential Vn{x) satisfies 
a higher order stationary g-KdV, g = 2n + 1, (Novikov) equation being a nonlinear ordinary 
differential equation with a linear highest derivative d^Vn/dx^ term [341 135j . (j2.ip belongs then to 
a class of finite-gap, or algebro-geometric system^. 

In a generic case, the emergent finite-gap potential Vn{x) has, however, singularities on M. An 
appropriate choice of the free particle non-physical eigenfunctions, 

, \ cosh Kjix + Tj), j = odd „ 
[smhKj(x -|- Tjj, J = even 

of energies = —nj, j = 1, . . . , n, gives rise to a nodeless Wronskian Wn{x). A non-singular 2n- 
parametric potential Vn = Vn{x; ki, ti, . . . , Kn,Tn) corresponds then to a reflectionless (Bargmann) 
system i7„ with n + 1 singlet states, separated by n gaps, n of which, of energies Ej^'^ = —Kp 
j = 1, . . . ,n, are the bound states, while the singlet state of zero energy, E = 0, lies at the bottom 
of the doubly degenerate continuous spectrum with E > 0. From another perspective, reflectionless 
potential Vn{x; ki, ri, . . . , Kn, Tn) describes n-soliton solutions of the KdV equation. 

Eigenstates of Hn'4>[n\ A] = Xil)[n; A], different from the physical bound states, are generated 
from eigenfunctions A] of the free particle, HqiI^[Q; A] = AV'iO; A], A / 

W(V>i,...,V>^,^[0;A]) 

where V'j are given by Eq. (j2.3p . Physical singlet bound states of Hn with A = — k| are obtained 
by the same prescription (12. 4p under the choice ViO; A] = sinhKj(a:; -|- Tj) for odd j, and 'i/'[0; A] = 
coshKj(a:: -|- Tj) for even j. The lowest singlet state of the continuous part of the spectrum of Hn 
corresponds to the eigenstate V'P; 0] = 1 of Hq. 

Transmission scattering amplitudes a[n;A;] for the continuous part of the spectrum E = fc^, 
/c > 0, of reflectionless system are defined by the scaling parameters k-,- [4], 

j=i -I 

The states (j2.4p have an alternative but equivalent representation, ip[n\X\ = j4„ . . . ^i^/^[0; A], 
generated by an n-sequence of the first order Darboux transformations, 

V-ijiA] =V'[('^,t)(,-);A] =^,-V'[j-l;A], (2.6) 

where («;, r)(j) denotes the set of 2j parameters ki, ri, . . . , Kj, r^-, and Aj = Aj[{K, are the first 
order differential operators defined recursively in terms of the states (j2.3p by 

did 



^Finite-gap periodic systems are given by the Its-Matveev representation of the form (|2.ip but with W{x) sub- 
stituted for Riemann's theta function [36]. If such a real periodic potential is regular on R, the Schrodinger (Hill) 
operator has a finite number of gaps in the spectrum. (|2.ip corresponds to the infinite period limit of the periodic or 
almost periodic finite-gap systems. An n-gap system also is characterized by the existence of the differential operator 
of order 2n + 1, related with a higher order Novikov equation, commuting with a Hamiltonian, see below. For the 
theory of finite-gap systems including historical aspects, see [121 137| . 
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The first order operator Aj annihilates the state Aj^i . . . Aitpj, that is a nonphysical eigenstate of 
Hj^i of eigenvalue —k'j. As inverse to ()2.6p . there is, up to an overall multiplicative constant, a 
relation 

ij[j-l;X]=A]^p[j;\]. (2.9) 

The zero mode of the first order operator A^ is l/{Aj^i . . . Aii/jj). It is the ground state of Hj of 
the energy —^j- 

A reflectionless j-gap Hamiltonian Hj admits two factorization representations 

Hj = ^U^j+i - 4+1 = AjA] - k] . (2.10) 

In particular, the free particle 0-gap Hamiltonian Hq = —-^ has an alternative representation 
Hq = ^1^1 — k\. From ()2.10p there follow intertwining relations 

AjHj_i = HjAj, A]Hj = H,^iA], j = l,...,n. (2.11) 

Let us take now a pair of n-gap reflectionless systems, 

Hn = Hn{Ki,Ti,...,Kn,Tn) and H'^ = HniK,'i, t[, . . . , K^, T^) ■ (2.12) 

Two sets of parameters are supposed to be completely different, or may partially coincide. If 
the two sets of the scaling parameters Kj, j = 1, . . . ,n, and k'-,, j' = 1, . . . ,n, do not coincide, 
the two systems have not only different spectra of bound states, but in accordance with ()2.5p . 
their transmission amplitudes are also different. If, moreover, nj ^ k!-, for all = l,...,n, 
all the energy levels of bound states for two n-gap reflectionless systems are different, and their 
transmission amplitudes are given by rational functions of k with different zeroes and poles. Having 
in mind that the factorization relations (I2.10p and the associated intertwining relations (j2.1ip are 
reformulated in terms of supersymmetric quantum mechanics construction, one can put a question: 

• What a supersymmetric structure is associated with reflectionless pair (I2.12p in a completely 
non-isospectral caseH characterized by inequalities kj ^ k'-, for all j, j' = 1, . . . ,n? 

Such a kind of supersymmetry of finite-gap pairs of refiectionless systems was not investigated 
yet in the literature. Instead, supersymmetry of the pairs {Hj, Hj^i), I > 1, belonging to the 
same Darboux chain (|2.1ip is usually considered. In particular, the pairs of reflectionless Poschl- 
Teller systems, see below, appear in the context of shape-invariance [381 EH E], they also emerge 
in the infinite-period limit of finite-gap crystal structures |20^ [22] . Supersymmetry of refiectionless 
Poschl-Teller pairs {Hj, Hj^i) was studied recently from the perspective of AdS/CFT holography 
and Aharonov-Bohm effect |40j . 

A special choice of the parameteres 

K'j = l^'j = jX', Tj=T, Tj=T', j = l,...,n, (2.13) 

results in two copies of the Ji-soliton potentials Vn = —n{n + l)K'^sech.^K{x + r) and = —n{n + 
\)K?seQh^ k{x + r'), which describe two mutually shifted refiectionless Poschl-Teller systems with 
n bound states. Since the partner potentials under the choice (I2.13P have exactly the same form, 
this corresponds to a particular case of a shape-invariance, whose analog in the case of periodic 
supersymmetric systems was called by Dunne and Feinberg ^self-isospectrality^ [15]. The exotic 

^ Using this term we neglect the fact that the continuous (scattering) parts of the spectra of the partner systems 
are the same, E > 0. 
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nonlinear supersymmetry of the simplest isospectral pair {Hi, H[) with ki = k'^, ti ^ r( was 
investigated and applied for the description of the kink and kink-anti-kink solutions of the Gross- 
Neveu model |41l [22] . One can expect that the self-isospectral pair of n-gap reflectionless Poschl- 
Teller systems with n > 1 should also be described by some not studied yet exotic nonlinear 
supersymmetric structure. 

In a more general case of the choice Hj = k'j, j = 1, . . . ,n, different from (12.13p . the partners with 
T 7^ r', r = (ti, . . . ,t„), are completely isospectral, their bound states energies and transmission 
amplitudes coincide, but the potentials have different form. We then arrive at the natural questions 
related to that formulated above: 

• How the supersymmetric structure of a general, non-isospectral case detects the coincidence 
of some of the scaling parameters of two systems in (j2.12p ? 

• Particularly, for a partial coincidence of the bound states energy levels, does the supersym- 
metry distinguish the coincidence of the scaling parameters of the same level, kj = k'j, from 
that corresponding to the case when distinct levels, Kj = k!-, with j ^ j' , coincide? 

• Is the case of a complete isospectrality of the two systems, kj = k'j, j = 1, . . . ,n, detected 
somehow by supersymmetric structure? 

• Does the case of self-isospectrality possess some special characteristics from the viewpoint of 
supersymmetry in comparison with a general case of isospectral systems with different form 
of potentials? 

In what follows, we study a peculiar supersymmetric structure of the pair (|2.12p . and, particu- 
larly, respond the highlighted questions. 

3 Supersymmetry of one-gap reflectionless pair with distinct scal- 
ings 

We first investigate the supersymmetric structure of the extended one-gap system 



described by the pair of reflectionless Poschl- Teller Hamiltonians Hi = Hi{k, t) and H[ = Hi{k,', t') 
with K ^ k' and arbitrary displacement parameters r and r'. This will allow us to trace how the 
restructuring of supersymmetry happens in the self-isospectral and to form a base for 

further analysis for n > 1, where we will restore index 1, omitted here to simplify notations, in the 
scaling and translation parameters. 

The choice of a non-physical eigenstate iPi{k,t) = coshK(2; + r), k > 0, r € M, of Hq produces 
a Hamiltonian of one-gap reflectionless Poschl- Teller system 



and first order operators Ai and A[ defined by Eq. (j2.7p . Operators Ai and A[ factorize the shifted 
for an additive constant Hamiltonians Hq and Hi, 




(3.1) 



Hi 



(P 2^2 



(3.2) 



dx^ cosh? k{x + t) 



Hi = AiA 




Hn = A\A 




(3.3) 
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and intertwine them, 

A\Hi = HoA{ , AiHq = HiAi . (3.4) 

A doublet of eigenstates in the continuous part, E = k'^, /c > 0, of the spectrum of Hi is 
constructed from the free particle plane wave states, 

= ^i(K,r)e±*'=^ = {±ik - Ktanh k{x + T))e^''''' . (3.5) 

The lowest singlet state with E = corresponds to a boundary case /c = of p.Sp . 

V'5' = tanhK(x + r). (3.6) 

Another, bound singlet state 

2 

ipi'^ =KsechK{x + T) (3.7) 
of energy E = —k^ is obtained from the doublet partner, 'tpi[K,,T) = sinhK(2; + r), of non-physical 

2 ~ 

eigenstate ^i(K,r) = coshK(2; + r) of Hq, ip^'^ {k,t) = j4i(k, r)^i(K, r). 

Based on intertwining relations (13. 4p and their analog for the system H[ = Hi{k',t'), we 
construct the second order operator 

Y2 = Y2{k,t,k',t')=Ai{k,t)A\{k',t') = AiA'^\ = Y2{k',t',k,t)=Y^, (3.8) 

that intertwines the partner Hamiltonians of the extended one-gap system (jS.ip . Y2H[ = H1Y2. 
Taking into account that Hq has an integral p = —i-^, one can obtain yet another, third order 
intertwining operator, 

Xs = Xs{k,t,k',t') = Ai-^A'i\ xI{k,t,k',t') = -X^i^' ,t' , k,t) = -X^ , (3.9) 

X3H[ = H1X3, which is independent from the second order intertwiner Y2. 

Intertwining relations in the reverse direction are obtained by a change t ,t' , that 

corresponds to a Hermitian conjugation of the corresponding relations, Y2H1 = H[Y2 , X-^Hi = 
H[xl, see Fig. HJi. 

The free particle integral P = —i-^ and intertwining relations (j3.4p also generate a nontrivial 
integral for the one-gap Poschl- Teller subsystem Hi{k,t), 

Z3 = Z3(k, t) = M-^A\ , Zl = -Z3 , (3.10) 

and the analogous integral, Z'^ = A[-^A'^, for Hi{k',t'). Integral ()3.10p is a nontrivial operator of 

a Lax pair for stationary KdV equation in the non-periodic case. 

Here and in what follows, the odd and even order intertwining operators are denoted by X 

and y, respectively, while the odd order integrals of the corresponding reflectionless systems are 

denoted by Z; the lower index indicates the differential order of these operators. 

Integral ()3.10p detects both physical singlet states of Hi {k, t) by annihilating them Z^ip^ {k,t) = 
2 ~ 2 

Z^'ip^'^ (k, t) = 0. The third state of its kernel is a non-physical eigenstate Vf*^ i-^) ~ 

V'l"'' (x)/dx/(Vr'' {x)f of Hi of energy —k , which is a linear combination of the physical bound 

2 

state Tp^"^ (x) of the same energy and of a non-physical eigenstate 'ipi{K, r) = coshK(x + r) of Hq. 

The extended system (13. ip has an obvious integral of motion 173. The intertwining relations 
together with integral (jS.lOp allow us to identify the nontrivial Hermitian integrals for the system 
-Hi, 

Ql;l = (^^t ' 2l;2 = ^f^sQl;!, Si-1 = , Si-2 = ^O-S^l;!, (3.11) 
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a) b) 

Figure 1: a) Non-isospectral one-gap Hamiltonians Hi (blue dot) and H[ (white dot) are intertwined 
by the second, Y2 and Y2 , and the third, X3 and X^, order Crum-Darboux operators via a virtual 
translation-invariant free particle system Hq (half blue/half white dot), b) In the isospectral case 
direct 'tunneling' channel for intertwining by the first order operators Jti and is 
opened. In both cases, Lax integrals Z3 and Z^, being the dressed forms of the free particle integral 
are the 'self-intertwining' generators for Hi and H[. 



= ( |/ ) > ^i;2 = '^s^i;! • (3.12) 

As erf = 1, we can take the integral F = (T3 as a Z2-grading operator. It classifies then Vi-a^ 
a = 1,2, as bosonic integrals, [oz-iT^i-.a] = 0, while the integrals (j3.1ip are identified as fermionic 
supercharges, {(T3, Qi-a} = {ca,^!;^} = 0, of the super symmetric structure of the extended system 
T-Li. There are other possibilities to choose F, which are based on reflection operators and classify 
the nontrivial integrals of the extended system in a way different from that prescribed by the choice 
F = (73. The alternative choices for F find some interesting physical applications, see [20 ^ l22 ^ HT | l42j. 
and we return to the discussion of this point in the last Section. 

Operators (|3.1ip and (13. 12^ are the Darboux-dressed integrals of the extended system described 
by the Hamiltonian Hq = d\a,g[HQ, Hq) composed from two copies of the free particle Hamiltonian 
Hq. The system Hq possesses the set of 2 x 2 matrix Hermitian integrals 

Io = cra, eabCTbP, Ip, cr-ip, a = 1,2. (3.13) 

The Darboux dressing, 

Zi = ViIqV\, Pi =diag(Ai(A€,r),Ai(K',T')), (3.14) 

transforms them into the integrals (j3.1ip and (|3.12p of Hi. 

We find the super algebraic structure of the system Hi by employing the intertwining and 
factorization relations (j3.4p and (j3.3p . It is given by the following nontrivial (anti)-commutation 
relations: 

{Qa.Qb} = 26abri{Hl,K)Fi{Hl,K'), {Sa,Sb} = 26abHlFi{Hl, k)Fi{Hi, k') , (3.15) 

{Sa, Qb} = 2eabFi{Hi,K)Vi , (3.16) 

[VuSa]=iHlF^{Hl,K,K')Qa, [VuQa] = -iF^{Hl,ti,K')Sa, (3.17) 

[V2,Sa]=iHlF+{Hl,K,K')Qa, [^2, Qa] = "^^ (^1, ^O^a , (3.18) 
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where Pi(?^i,k) ='Hi + k^-1, Pi(?^i,/C) =Hi + lO,lC = diag(K',K), 

Pq (-Hi, k') = PiCHi, k) - Pi (7^1, k') = (k^ - . 1 , (3.19) 

P^-CHi,^,^') = PiCHi,^) + Wi{ni,K') = 2Hi + (k^ + k'2) • 1, and to simplify the formulae, we 
omitted the index n = 1 in the supercharges and bosonic integrals. Though in the final expression 
for Pq in (j3.19p the dependence on l-Li disappears, it is indicated here in the arguments having 
in mind a further generalization for the n-gap case, where this structure is substituted for the 
polynomial of order n — 1 in Hamiltonian. 

The one-gap extended reflectionless system (13. ip is described therefore by a nonlinear super- 
algebra generated by four fermionic supercharges, Qi-a and 5i;a, and by two bosonic integrals H, 
■Pi;a. The fermionic integrals are constructed from the intertwining operators of the second and 
third orders, whose composition produces nontrivial third order integrals of Lax pairs of the one 
gap non-isospectral subsystems. In this super symmetric structure, Hamiltonian plays a role of the 
multiplicative central charge. The nonlinear superalgebra depends here on the scaling parameters 
K, and k! via the polynomials Pi , and Pq , but does not depend on the displacement parameters 
r and r'. 



4 Supersymmetry of a one-gap self-isospectral pair 

For the isospectral extended system T-ii with k = k', the partner potentials have the same form 
and are mutually displaced. This one-gap self-isospectral case is special from the viewpoint of 
supersymmetric structure. As follows from ()3.17p and (I3.19p . for k = k' the integral 'Pi;i, composed 
from the third order integrals of Lax pairs of superpartner subsystems, commutes with all the 
integrals, and so, transmutes into a bosonic central charge of the nonlinear superalgebra. We show 
now that the supersymmetric structure in this case undergoes even more radical changes. 
For K = k' the following reduction takes place: 



XsiK, T, K, t') = {Hi{k, t) + Xi{k, T, t') - C{k, T - T')Y2iK, T, K, t') , (4.1) 



where 



Xi(k,t,t') = utanh k{x + t) + Ktanh k{x + t') + C{k.,t — t') (4.2) 

dx 

= Ai{k,t)-A{{k,t')+aI{k,t-t'), (4.3) 

Ac{k,T - t') = +C(k,T - t') , C(k, T - t') = KCOthK(T - t') . (4.4) 

dx 

Relation (14. ID means that for r ^ r', the first order operator Xi = Xi{k,t,t') should be taken 
as a basic odd order intertwining operator instead of X^{ 

XiHi{K,T') = Hi{K,T)Xi, xI{k,t,t') = -Xi(k,t',t) = -X[. (4.5) 

Note that in the limit r' — ?■ ±00, we have H[ Hq and Xi — ?> Ai, while for r ±00, Hi — )• Hq 
and Xi — 7> —A'l- This is coherent with the intertwining relations (|3.4p . 

Because of (14. ip . the third order integrals Si-a are reducible, Si-a = {Hi + ^<^)Si-a — CQi-a , and 
have to be changed for the first order irreducible integrals 

■51;! = f ' '^i;2 = io-aSi-i ■ (4.6) 



''There are four bosonic integrals if one counts the integrals "Hi and 0-3. 
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Integrals Si-a correspond, in accordance with (j3.14p . to the dressed form of the integrals 
Sa = ^abO'bP + Caa, of the extended free particle system T-Lq = diag(ffoi Hq), VsaV^ = Si-a{T-Li + k^). 
Alternatively, the first order matrix operator si = + Cai, or S2 = ic^si, can be considered as a 
first order Hamiltonian of the free Dirac particle of mass \C\ in (1+1) dimensions, while its dressed 
form, Si-i, can be identified as a Bogoliubov-de Gennes Hamiltonian describing the kink-antikink so- 
lution in the Gross-Neveu model [30]. Function A(^, X) = k (tanh (.^ — A) — tanh + A) + coth 2A), 
that appears in the structure of Xi with ^ = k(x + ■^-^—) and A = —k '^~2 > then a sense of a 
gap function [21]. 

The following relations are valid: 

X^xI = H^{k,t)+C\ (4.7) 

XiAi{k,t') = Ai{k,t)Ac{hi,t - t') , A\{k,t)Xi = Ac{^,T - t')A\{k,t') . (4.8) 
The employment of (|4.7p . (|4.8p together with (j4.5p gives nontrivial nonlinear superalgebraic relations 

{5l;a, Si,b} = 25abhc , {Ql-a, Ql;b} = ^S^bhl , (4.9) 

{Sl-a, Ql;b} = 25abCK + 2e^bPl;l , (4.10) 
[Vl,2, Si,a] = 2i{hcQl;a - ChJi-a) , [Pl-2, Ql-a] = 2iK{CQi,a " KSi,a) , (4.11) 

which substitute nontrivial superalgebraic relations ()3.15p . (j3.16p . ()3.17p and (jS.lSp of the general, 
non-isospectral case n = 1. Here we denoted /i^ = 'Hi + k'^, he = Tii+C^. As > k^, the spectrum 
of he is strictly positive, and the Lie sub-superalgebra generated by the first order supercharges Si-a 
corresponds to a broken N = 2 super symmetry. The Vi-i commutes now with all the supercharges 
in accordance with the observation made at the beginning of the Section. 

While the third order intertwining operator (j3.9p is well defined at k = k' , t = t' and reduces to 
the integral Z3{k, t) of Hi{k, t), the first order intertwining operator Xi{k, r, r') in the limit t' ^ t 
reduces to the operator ^ shifted for an infinite additive constant term ±00 in dependence on which 
side the difference (r — r') tends to zero. In this case extended Hamiltonian (j3.ip reduces just to the 
two identical copies of the Poschl- Teller Hamiltonians, T-Li{k,t) = diag{Hi{K,T), Hi(k,t)). The 
integrals Si-a, (2 = 1)2, can be renormalized multiplying them by 1/C{k, t — r'), and taking a limit 
t' — )• T. In such a way they are reduced to the trivial integrals ctq, a = 1, 2, of 'Hi{k, t). The second 
order intertwining operator (|3.8p reduces in the limit r' — )■ r to Hi{k, r) + k^, and the second order 
supercharges Qi-a are reduced to the same trivial integrals aa multiplied by a shifted for a constant 
Hamiltonian, Qi-a {'Hi{k,, t) + k? ■ l)o"a. The only nontrivial integrals we have in the limit t' ^ t 
are the bosonic third order integrals Vi-ai^jT). 

The special case of self-isospectrality in the one-gap extended system T-Li is detected, therefore, 
by a radical change of nonlinear supersymmetric structure. One of the bosonic integrals, Vi-i, turns 
into a central charge, and two third order supercharges are substituted for the supercharges of the 
first order. The reduction of the order of the half of the supercharges at k = k' originates from 
relation (|4.1j) and is accompanied by appearance of dependence of the superalgebraic structure 
on the distance between mutually shifted one-soliton partner potentials by means of a constant 
C = KcothK(r — t'). In other words, one can say that in a generic case k 7^ k', the Hi and H[ are 
intertwined by the third order operators and X^, side by side with the second order operators 
Y2 and Y2, via the free particle (zero gap) system, and the supersymmetric structure does not 
feel a relative distance r — r' between the corresponding one-soliton subsystems because of the 
translation invariance of Hq. For kind of a 'tunneling' channel is opened: the one-soliton 

subsystems are intertwined then directly by the first order operators Xi and xj, and the modified 
supersymmetric structure detects a 'tunneling distance' r — r', see Fig. [TJa. 
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5 Supersymmetry of an n-gap system: complete isospectrality 
breaking 

The discussion of supersymmetric structure for extended system composed from two subsystems 
having n > 2 bound states requires to distinguish three cases: 

• Complete isospectrality breaking, when Ki ^ k'j for all i,j = 1, . . . , n, with no restriction on 
displacement parameters Tj and rj. 

• Partial isospectrality breaking, in which some, but not all, scaling parameters Kj and k'j of 
the two subsystems coincide. 

• Exact isospectrality, that is characterized by the complete coincidence of the sets of the scaling 
parameters, k = k', accompanied by a restriction t ^ f' . 

The case of a complete isospectrality breaking for n > 1 is a direct generalization of that for 
n = 1 case with ki / which was studied in Section [3l It is discussed in the present Section. 
Other two cases are more involved. Though they generalize somehow the picture of the one-soliton 
case (n = 1) with investigated in the previous Section, the corresponding analysis 

for n > 1 requires a generalization of the described Crum-Darboux transformations scheme. New 
peculiarities appear there, and those two cases deserve a separate consideration. To understand 
the picture, we study the case of n = 2 in the next Section, and then in Section [7] the results will 
be extended for a generic case of n > 2. 

With these comments in mind, let us consider an extended system 



composed from a completely non-isospectral pair Hn = Hn{K,T) and H^-^ = Hn{K,' ,t') of the form 
(|2.12p . where k = {ki, . . . , Kn), r = (ri, . . . , r„), and we assume that there is no coincidence in the 
sets of the scaling parameters of the two subsystems, kj ^ k'-, for all j, j' = 1, . . . ,n, see Fig. [2^. 

Following the general picture described in Section [2l Hamiltonian Hn = Hn{K,f) can be inter- 
twined with a free particle Hamiltonian Hq by order n differential operators A^^ = An(K,'T) and 




(5.1) 





(5.3) 



Making use of these relations, we construct an order 2n operator 




/ 1 




(5.4) 



■n 



where A' 



■n 



A„(K',r'), and two operators of the order 2n + 1, X2n+i and Z2n+i 



X2n+l(K,f;K',f') = A„— a; 



■n 



2n+l ~ 



X2n+l{K',T';K,T) = -X'2 



2n+l ) 



(5.5) 




'2n+l ~ 



(5.6) 
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Figure 2: a) An n-gap pair with complete isospectrality breaking. Each subsystem, Hn and H'^, 
is specified by indicating the set of intermediate, virtual, systems in the plane k — t via which the 
edge points are connected to the free particle by means of the first order Darboux generators Aj 
and j4j, not shown here. Figs, b) and c) illustrate two alternative representations for the same 
n = 2 system, that is related to the two different factorizations of the second order Crum-Darboux 
generator A2. In the case b) the virtual system is regular, while in the case c) it is singular. So, a 
system is specified not only by indication of the set of points in the k — t plane, but also by the 
path via these points to a free system Hq. 




Operators Y2n and X2n+i intertwine the components of the matrix Hamiltonian Tin □, 

y^nH'n = HnY2n , X2n+lH'n = HnX2n+l , (5.7) 

while Z2n+i(K, r) is an integral for r), 

[Z2n+l,Hn]=0. (5.8) 

Taking into account that the coefficients of the (2n + 1) order differential operator Z2n+i may be 
expressed in terms of the potential and its derivatives of the order less than 2n + 1 [35] , relation 
()5.8p means that the potential Vn satisfies a higher stationary g-KdV equation with g = 2n + 1, 
mentioned in Section [2j 

In correspondence with an identity ^In+i — ~Hn 11^=1 (-^n + ^j)^i integral .^2n+i detects all 
the physical singlet states of Hn of energies E = Q and Ej = —Kj by annihilating them. These are 
constructed from the free particle singlet eigenstate ipQ = 1, = and non-physical partners 

of the states (12. 3p . Vi = sinhKi(x + ri), ip2 = coshK2(x + T2), . . . , ipn ' = ^n^Pj, j = 1, • • • ,n. 

Other n states of the kernel of .^2n+i sue non-physical doublet partners of the bound states ipn ^ , 

V'n '{x) ='4)n '{x) f dx/{lpn ' (x))^. 

With the described operators, we construct six matrix integrals Qn;a, Sn-a and Vn-a for the 
extended system Tin in the form similar to that in (j3.1ip and (j3.12p by changing Y2, X3 and for, 
respectively, Y2n, -^^271+1 and .^2n+i- As in the n = 1 case, these integrals correspond to a dressed 
form of the integrals of the extended free particle system T-Lq obtained by means of Eq. (13.14p with 
the change of Pi for Vn = diag (A(k, f), A„(k', r ')). 

Applying factorization and intertwining relations, and products of corresponding generators 
collected in Appendix, we find that the superalgebra ()3.15p . (I3.16p . ()3.17p . (I3.18P of the one-gap 

^Intertwining relations through multi-step ladders of linear Darboux generators and their superalgebraic reducibil- 
ity have been recently reviewed in [43], but in a very general and abstract form. 
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case is generalized for 

{Qn;a, Qn;b} = 25„f,P„(?^„„ k)P„(?^„, k') , {Sn-a, = 25^b?^„P„(?^„, k)P„(?^„, k') , (5.9) 

{5„;a, Qn;b} = 2e,feP„ (?^„ , £)P„;l , (5.10) 
{Vn;l,Sn-a]=innK^l{nn,K,K')Qn;a, [Pn;!, Qn;a] = -iF-^iCHn, K, K')Sn;a , (5.11) 

[Vn;2,Sn;a] = lUnKi^n, K, K')Qn;a , [K;2, Qn;a] = -iP+(^n, ^')'5n;a , (5.12) 

where P„(?^„, ^) = n"=i('Wn + ■ 1), P+(^n, ^, ^') = F„CH„, ^) + P„(?^„, k'), P;:_i(-Hn, K, Ac') = 

Fni-Hn, fi) - Pn(^n, k'), Pn(^n, /C) = ni=l(^n + /C|), /Cj = diag {k'j, Kj). 

Operator F~_i{'Hn, k') is a polynomial of order n — 1 in the extended Hamiltonian Hn that 
vanishes for k = k' . Then Eq. (|5.1ip signals that the super symmetric structure of the n-gap 
reflectionless system Hn with exact isospectrality simplifies as in the case n = 1: the integral 
Vn,! turns into bosonic central charge of the nonlinear superalgebra. Moreover, from the form of 
polynomial in "Hn coefficients in superalgebra, one can expect that the super symmetric structure 
should undergo some radical changes even in the case when not all the pairs of the scaling parameters 
coincide but only part of them. For instance, if k'^, = Kj for some indexes j' and j, which may 
coincide, j' = j, or may be different, j' ^ j, the same factor {T-ln + ■ 1), or its square, appears 
in all the structure coefficients of the superalgebra. By analogy with the n = 1 case this indicates 
that some fermionic supercharges may be substituted for supercharges of a lower differential order. 
To understand what changes the supersummetric structure undergoes in the cases of a partially 
broken or exact isospectrality, we investigate in detail the extended system ()5.ip for the case of 
n = 2 in the next Section. 

6 Supersymmetry of the n = 2 extended system 

Explicit form of the supersymmetric structure for extended n = 2 system with completely broken 
isospectrality follows as a particular case from a generic consideration of the previous Section. 
Before analyzing the partially broken and exact isospectrality cases, we first discuss some properties 
of the n = 2 reflectionless system of the most general form. It is a particular case of such a system, 
described by the two-soliton Poschl- Teller Hamiltonian, that appears in ^p^ field theoretical model 
with a double well potential, where it controls the stability of the kink and anti-kink solutions. 

6.1 Generic ri = 2 reflectionless system 

Explicit form of the Hamiltonian of two-gap reflectionless system of a general form is 

H2{k,t) = -— + V2{x-k,t), (6.1) 
V2{x;k,t) = — 2(^2 — )~^ (ft;2csch^«;2(a; + T2) + Kisech^Ki(rE + Ti)) tt;^(x; K, r) , (6.2) 

where 

w{x] K, r) = {k1 — K2) {k2 coth K2{x + T2) — Ki tanh ki{x + ti))~^ . (6.3) 
In the limit T2 ±00, the two-soliton system (j6.ip transforms into that of the one-soliton case, 

V2 -2Klsech^Ki{x + ri T 6) , (6.4) 

where a shift parameter is defined by a relation sinhKi^i = ki/ K2 — nf. In another limit, ri — )• 
±00, the two-soliton potential transforms into the one-soliton potential given by an expression of 
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the form (j6.4p but with the index 1 in the parameters changed for 2; the shift parameter ^2 is given 
then by a relation sinhK2.^2 = ^2/\/^2 ~ '^i- '^^^ indicated hmits correspond to a picture of a 
two-sohton scattering described by the KdV equation, where the 'one-gap soHtons' of amphtudes 
2k1 and 2^^ in such a process suffer asymptotically only temporal shifts [2]. 

Singlet bound states of the system (|6.ip . V2 ^ — j = 1)2, of energies E = —k\ and 

E = —K2 are obtained from the doublet partners, ipi = sinhKi(a; + ri) and ^2 = cosh ^2(2; + T2), of 
non-physical eigenstates ipi = coshKi(x + ri) and ip2 = sinhK;2(x -|- T2) of Hq by applying to them 
the second order composite operator 

A2{fi,T) = A2{K,f)Ai{ni,T) , (6.5) 



2 2 

^2 ^ = KisechKi{x + Ti) w^x; k,t) , ^^2 ^ = — «^2csch ^2(2; + T2) ^^(x; k, r) . (6.6) 
A2{k,t) = {Ai4'2) -^Yir~r\ = -^i(«^i)n) + t/^(2;;K,f) , (6.7) 



Here 



and Ai is defined by Eq. (|2.7p . Function (j6.3p satisfies the identities 

du'/(ix = iV'2, (6.8) 
w"^ + 2Kiu;tanhKi(a; -|- ri) = ^V2 + k.^ — kI , (6.9) 
w'^ + 2k2W coth K2ix + r2) = ^V2 + k{- kI, (6.10) 

which will play a fundamental role in what follows. 

The doublets of the states of the continuous part of the spectrum with E = k'^ >Q are obtained 
from the plane wave states of the free particle, 'ip^'^ = A2e^*'^^, 

^f^= [-(fc2 + K?) + (±iA;-KitanhKi(2; + ri))u;(x;K,f)]e^*'=^. (6.11) 

The boundary case A: = gives a singlet, zero energy edge state 1/^2 the bottom of the continuous 
spectrum. 

The particular case of refiectionless n = 2 Poschl- Teller system, 

H2{k, t) = — —J — 6K;^sech^K(x -|- r), 

is obtained by putting K2 = 2ki = 2k and T2 = n = r. In this case, the function (16. Sp and the 
operator (|6.7p are reduced io w = — S/ttanhx and ^2 = ^ — 2Ktanhx, the indicated bound states 
are transformed, modulo overall multiplicative constants, into '^2''^ = sinh^sech^X (-^ = ~^^) ^'^d 
1^2'^'^ — sech^x (-^ = ~4k^), while the zero energy singlet state is ■02 = ^ ~ Stanh^X) where we 
use the notation x = i^{x + t). 



6.2 Generalized Crum-Darboux transformations scheme 

We have constructed a generic n = 2 refiectionless Hamiltonian (16. ip by employing a sequence of 
two Darboux transformations described in Section [21 namely, by using first the non-physical free 
particle state tpi = coshKi(x + ti), and then the state 11)2 = sinhK2(x + T2). The same final result 
also can be achieved with the interchanged order of the indicated states. This corresponds to the 
alternative factorization of the second order operator (16. 5p . 

A2 = S2S1 , (6.12) 
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which intertwines H2 with the free particle Hamiltonian, A2i/o = H2A2, A2H2 = HqA2, see Fig. 
[2]3, c. The first order operators Bi and B2 are obtained from Ai and A2 via the substitution 

IT vr 

Kl^K2, Ti^T2 + t- =T2, T2 Ti + I- = Ti . (6.13) 

2k,2 2ki 



This substitution leaves invariant the Hamiltonian (j6.ip . the second order intertwining operator A2, 
and the function (|6.3p . It also leaves invariant the states (j6.1ip of the continuous spectrum, including 

the singlet edge state of zero energy, but interchanges the bound states (j6.6p . V2 ^ ~^ "^^2 ^' 

ip2 ^ ^ 'i'4'2 ^ ■ Transformation (|6.13p changes, however, the first order intertwining operators Ai 
and ^2, which are regular on M^, for the singular first order operators 

Bi = Bi{ti2,T2) = ^ ~ K2CothK2(2; +T2) , B2 = B2{k,t) = -i?|(K2,T2) +'w{x;k,t) . (6.14) 

In terms of the first order operators (j6.14p we have = Kg, Hi = Hi{k2, f2) = B\B\ — \^ = 

bIB2 - kI, and H2 = B2BI - nj. This means that with the alternative factorization (|6.12p . the 
operator A2 intertwines Hamiltonian (j6.ip with Hq via the n = 1 system described by a singular 
Hamiltonian ^ ^ 

H, = Hi{K2,f2) = + ■ ^2 ^? , ^ • (6-15) 

OX^ Smh K2(X + T2) 

In what follows, singular Hamiltonian ()6.15p will appear only as a virtual, or intermediate system, 
and the described generalization of the Crum-Darboux scheme will allow us to identify nontrivial 
intertwining operators for n = 2 extended system with partially broken and exact isospectrality. 
The picture with the alternative factorizations generalizes for the case n > 2. In this context 
it is worth to note that the change of the order of the free particle non-physical states (|2.3p in 
the construction of a reflectionless system H^, in comparison with that described in Section [21 
corresponds to a certain permutation of the columns of the Wronskian ()2.2p . This produces no 
effect for potential in equation (j2.ip . 

To conclude the discussion of the generalized Crum-Darboux transformations scheme, we present 
here the relations which are helpful for computation of the corresponding superalgebraic structures, 

Xi{K,f,f')B[ = BiAc{K,T-T'), BlXi{K,f,f') = Ac{K,T-T')B[\ (6.16) 
Xi {k, f , r') A[ = B^AciK, t-t'), B^X^ {k, f, r') = Ac{k, f - t')A'^ , (6.17) 
Ac{ti,T-f')B'^ = A\Xi{K,T,f'), (6.18) 

where Ai = Ai{k,t), A[ = Ai{k,t'), Bi = Bi{k,,t), B[ = Bi{k,,t'), f = r + i-^, and f' = t' + 
These identities can be obtained from (|4.8p via the substitution (|6.13p . 

6.3 Generic case of partial isospectrality breaking 

Now we are in position to discuss the supersymmetric structure of the extended n = 2 systems with 
partially broken and exact isospectralities. We first consider three cases of partial isospectrality 
breaking, in which one discrete energy level — k| of the subsystem H2 coincides with any of the 
two discrete energy levels — of the partner Hamiltonian H2, but the corresponding translation 
parameters are different, Tj 7^ t'-,. All these cases are described by a similar supersymmetric 
structure. Then, in the next subsection, we analyze the superalgebraic structure of the same three 

cases but with coinciding associated translation parameters, r,- = r',. 

3 

We start with the case of partial isospectrality breaking characterized by the conditions 
• = k'i, Ti 7^ r(, K2 7^ ^2, no restrictions on T2, T2, (6.19) 
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a) b) c) 

Figure 3: The n = 2 pairs with partially broken isospectrality. 

see Fig. [3h. 

The subsystems H2 = H2{ki,ti, K2,T2) and H2 = H2{ki,t[, k^^t!^) of the extended matrix 
Hamiltonian 'H2 are related by irreducible intertwining operators of orders 4 and 3, 14-^2 ~ -^^2^; 
yIH2 = H2YI, X^H2 = H2X^, X^'^H2 = H2X^'^. Y4 is given, in correspondence with the generic 
form (I53D, by Y4, = A2A'^ , while 

X^ = A2{Ki,K2-,n,T2)Xi{Kl-,n,T[)Al{Ki,K'2-,T[,T^) = A2XiA'2^ (6.20) 

appears instead of the fifth order intertwining operator X^ = A2^A2 because of the reduction 

X, = {H2 + Kl)Xi - C{Ki,Ti - . (6.21) 

As follows from (j6.20p . the reduction (j6.2ip is related to the opening of a 'tunneling channel' via 
the virtual isospectral pair of one-gap systems //i(ki,ti) and Hi{ki,t[). 

Taking the products of the described intertwining operators and Lax integrals of order 5, 
^5 = A2^A^, Z5 = A'2^A2^ [^5,i^2] = 0, [^5,-^2] = 0> presented in Appendix, we find the 
super algebraic structure of the system 7^2 with partially broken isospectrality (|6.19|) . It is dis- 
played below in the form that unifies (|6.19p with two other similar cases. 

A partial isospectrality breaking with coinciding ground state energy levels, 

• K2 = 1^21 ■^a 7^ 'T21 1^1 A i^'ii rio restrictions on ti, t[ , (6.22) 

is similar to the previous case, see Fig. [3)3. Intertwining operator I4 and integrals Z5 and Z'^ are 
given by generic formulae with restriction (I6.22p . The third order irreducible intertwining operator 
can be obtained from (|6.2Up via the substitution (|6.13p . 

Xi = B2Xi{K2;f2,f^)B'2^ , (6.23) 
where B2 and B2 are given by Eq. (|6.14p with K2 = K2, while 

^1(1^2', t'2,t'2) = ^ - COthK2(2; + r2) -I- K2 COthK2(2; + T"2) + C(^2,''"2 - T2) 

= BiiK2,T2)-Bl{K2,T^)+AliK2,T2-r^). (6.24) 

Though all the three first order operators that appear in factorization of X3 in ()6.23p are singular, 
the third order intertwining operator itself is regular on M}. This follows just from the reduction 
relation for the fifth order intertwining operator for the case ()6.22p under consideration, 

X5 = {H2 + Kl)Xi - C{k2,T2 - T^)Yi . (6.25) 
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The third order intertwining operator (j6.23p reahzes the intertwining between H2 and H2 by means 
of a 'tunnehng channel' via a pair of singular one-gap Hamiltonians Hi{K2,f2) and Hi{K2,f2) of 
the form (l6T3|) E|. 

The supersymmetric structure for partial isospectrality breaking 

• ^vi = K2, K2 7^ k'i, no restrictions on ti^2, t[ 2 1 (6.26) 

see Fig. [3j:, is generated in a similar way. Here, the third order irreducible intertwining operator is 

Xi"" = A2Xi{Ki;Ti,f!,)B'2^ , (6.27) 

where B2 = B2{k[, K2,t[,T2) is given by Eq. (j6.14p . and f2 = T2 + ij^- In this case, we have a 
reduction relation 

X5 = {H2 + nDxi"" - C(ki, n - f'2)Yi . (6.28) 

Unlike the two previous cases, C(ki,ti — fg) = KitanhKi(Ti — Tg) is regular for any values of ri 
and T2 associated with coinciding scaling parameters [§. 

The superalgebra for the described three cases of partial isospectrality breaking can be presented 
in a unified form 

{Sa: Sb} = 25abhdhd'hQ , {Qa, Qb} = "^Sabhf hdhd' , (6.29) 
{Sa, Qb} = 25abCihihdhd' + 2eabhd',dVi , (6.30) 

\Vu Sa] = - ^'!){hc,Qa - ClhSa) , [Vl,Qa] = i{4 - K'j)h^{ClQa " hSa), (6.31) 

[V2, Sa] = i{hd + hd'){hciQa - CihiSa) , [V2, Qa] = ^ih-d + hd')hi{CiQa - hiSa) . (6.32) 

Here hi = T^a + ^f, hd = H2 + i^% hd' = ^^2 + ^^ , hd',d = ^2+diag(Kf , k^), = 'H2+Cf, I = 1,2,3, 
Ki denotes the coinciding scaling parameter of the pair, and k'^ correspond to other, not coinciding 
scaling parameters of the subsystems H2 and H2, respectively, while Ci = C{ki,ti — t[) for (j6.19p . 
C2 = C{k2,T2 — T2) for ()6.22p . and C3 = C(ki,ti — f2) for the case (|6.26p . Notation Sa refiects the 
reduction S2-a = (^2 + >^^)Sa — CiQa of the supercharges constructed in terms of X^ and xl, and 
to simplify notations, we do not supply the supercharges with index /, and omitted the index n = 2 
in all the integrals. 

The fact of a partial isospectrality breaking is refiected here in the superalgebraic structure. 
On the one hand, relations (j6.29p . (j6.30p and (|6.32p are similar to superalgebraic structure (j4.9p . 
()4.10p and (14. lip of the n = 1 isospectral case. At the same time, the commutators in (I6.3ip . being 
of the nature of those in (j3.17p for the n = 1 non-isospectral family of the systems, show that a 
'non-centrality' character of the Lax matrix integral 1^2-1 is measured by the scale of isospectrality 
breaking, — k'?. 



6.4 Partial isospectrality breaking with coinciding associated translation pa- 
rameters 

Let us discuss now the supersymmetry of the systems with partial isospectrality breaking, in which 
one discrete energy level, kj = k'-,, and the associated translation parameters, Tj = rj,, coincide, 
see Fig. Hti, b, c. The two cases corresponding to either ki = k'^ or K2 = ^2 ^'^^ similar. For them, 

^By shifting the argument x ^ x + iS, where J is a real constant, one can translate all the consideration for the 
case of ■PT-symmetric quantum systems [35] with Hi and H[ to be regular isospectral Hamiltonians, see [46| . 

®The operator (|6.27p intertwines H'2 and H2 via the virtual one-gap systems H[ and Hi of different, singular and 
regular, nature. After the imaginary shift mentioned in the previous footnote, the latter pair will transform into 
regular one-gap Poschl- Teller PT-symmetric Hamiltonians. 
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b) 





Figure 4: The pairs with partiahy broken isospectrahty, in which the translation parameters as- 
sociated with the equal scahng parameters do coincide. In the case a), a common virtual system 
corresponds to a regular one-gap Poschl- Teller system. In the case b) such a common virtual system 
is singular. In the case c), the partners can be intertwined via a pair of n = 1 virtual systems, one 
of which is singular. 



supersymmetry undergoes restructuring, and is generated by intertwining operators of the second, 
Y2, arid fifth, X5, orders, and by the fifth order integrals and Z'^. The fifth order operators, X5 
and Z5, in this case include in their structure the third order integral of the corresponding common 
virtual n = 1 system. 

For the sake of definiteness, consider the case ki = k[, ti = r(, K2 7^ ^2- have = A2Z^A'2\ 
Z5 = A'^Z'iA'i^ ^ and Z'^ = ^2Z3^2^ where Z3 = Zs{ki,ti) = Ai{ki,ti)-^A\{ki,ti) is the third 
order Lax integral for the common Poschl- Teller virtual system Hi{ki,ti). The second order 
intertwining operator has a form = A2A2, with A2 = ^2(^1, K2,ti,T2) and A2 = {ki, Kg, ti, Tg), 
and the fourth order intertwining operator I4 = A2A2^ of a generic case reduces as 

n = {H2 + «?)y2^ . (6.33) 

The second order operator can be obtained also from the third order operator (I6.20p of the 
case ()6.19p considered above. Indeed, multiplying ()6.20p by —C~^{ki,ti — r{), and taking a limit 
r( Ti, we get Y^. So, the change of supersymmetric structure is related here to a singular nature 
of C{ki,ti — t[) in the limit r{ — > ri. Another case, with K2 = ^2, T2 = Tg, ki 7^ k'i is treated in a 
similar way, and super algebraic structure for these two cases can be presented in a unified form: 



{Sa, Sb} = 26ab'H2hihdhd', {Qa, Qb} = 26abhdhd', (6.34) 

{Sa, Qb} = 2eabhd',dVi , (6.35) 

[Pi, Sa] = - K'i)n2h^Qa , [Vl,Qa] = -^(^^1 " 4)Sa , (6.36) 

[r2,Sa] = iU2hi{hd + hd')Qa , [V2, Qa] = -i{hd + hd')Sa . (6.37) 



Notation Qa reflects here the reduction Q2-a = ('^2 + i<'i)Qa, and, again, we omitted the index 
n = 2 in the specification of nontrivial integrals. 

The case ki = K2, ti = T2, 1^1 1^2 different from the two previous ones because the 
corresponding parameter-dependent function C{ki,ti — f'2) = ki tanh ki(ti — Tg) is non-singular for 
any values of ri — T2, and, moreover, turns into zero at ti = Tg. Here, the intertwining operators 
are Y4, and given by Eq. (j6.27p with ri = T2. A non-singular nature of the latter is seen 

from (16.28p . The superalgebra for this case is obtained directly from (I6.29p - (l6.32p just by putting 
there C3 = 0. Though here the irreducible intertwining generators are different in comparison 
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with the previous two cases, the resulting superalgebra ()6.29p - (|6.32p with C3 = has a similar 
form to (j6.34p -( r6.37p . Notice also a remarkable similarity of (I6.34p - (]6.37p with the superalgebra 
(j3.15p - (|3.17p of the n = 1 non-isospectral case. 

We see that in all the three cases of partial breaking of isospectrality with corresponding coin- 
ciding translation parameters (associated with coinciding discrete energy levels), the superalgebraic 
structure does not depend on the two remaining translation parameters associated with the second, 
different discrete energy levels. 

In all the cases of partial isospectrality breaking described in this and previous subsections, the 
total order of the two basic intertwining operators is the same, 3 + 4 = 2 + 5 = 7, being less in 2 
in comparison with the complete isospectrality breaking case. 

6.5 Exact isospectrality with a common virtual n = 1 system 

The super symmetric structure of the systems with exact isospectrality, ki = k\ and K2 = ^2, 
depends on whether the corresponding translation parameters are different, Tj ^ rj, j = 1,2, 
or they coincide in one of the pairs The analysis of the second case, see Fig. [5ti, b, is more 
simple, and we first consider it supposing, for the sake of definiteness, that ri = t{, T2 7^ Tj. The 




a) b) c) 

Figure 5: The n = 2 isospectral pairs with a common regular, (a), or singular, (b), virtual system. 
A general case of the n = 2 isospectral pair with tj / rj, j = 1, 2, is illustrated by c). 

intertwining operators for such an isospectral system with a common regular virtual one-gap system 
Hi{ki,ti) are 

= , = B2X^ {K2,f2,r'2)B'^ , (6.38) 

where in A'2 and B'2 we assume that k'- = kj, j = 1,2, and t[ = ri, T2 ^ T2. They can be obtained 
here via the reduction relations of generic intertwining operators, 

X5 = {H2 + 4)Xi - C{k2,T2 - t'2)Y^ , = {H2 + nl)Y2^ . (6.39) 

The intertwining generators Y^ and X^, and the corresponding reduction relations for the exact 
isospectrality case Kj = k^-, j = 1,2, = T2, t[ ^ ti are obtained from (|6.38p and (|6.39p by 
changing ki f-)- K2) -H- T2, t{ -H- Tg, and A2 -v^ B2. 

The nontrivial relations of superalgebraic structure for the isospectral case with ti = r( , r2 7^ T2 

are 

{SaJb} = 26abhc,hl {Qa, Qb} = 2dabhl (6.40) 

''The case when both pairs of translation parameters coincide corresponds to 'H2 composed from the two copies 
of the same Hamiltonian H2 ■ Such a system 'H2 is described by a trivial supersymmetric structure to be similar to 
that discussed for 71 — 1 case in Section [l] with integral Zs changed for Z5. 
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{Sa, Qb} = 26abC2hih2 + 2eabVl, (6.41) 
[P2,Sa] = 2ihi{hc,hiQa - C2h2Sa), [V2, Qa] = 2i/l2 (Ca/ll Qa " /i25a), (6.42) 

where C2 = K2 cothK2(T2 — Tj), hi = T-L2 + nf, i = 1, 2, and hc^ = ^-2 + The superalgebra for 
the isospectral case with T2 = T2, ti t[ is obtained from the displayed one by changing C2 ^ Ci, 
hi -H- /i2 in the right hand side expressions. The supersymmetry (|6.40p . (j6.4ip . (|6.42p has the 
structure similar to that for the n = 1 isospectral case. 

As it is expected, the integral P2;i transmutes here into the bosonic central charge of nonlinear 
superalgebra. The total order of the basic irreducible intertwining generators reduces in two in 
comparison with the partially broken isospectrality case and equals the order 5 of Lax integrals Z5 
and Zg. In correspondence with this, the anticommutator of the second order, Q,2;a, and the third 
order, S2-a, supercharges taken with different values of indexes a and b is equal to the central charge 
1^2:1 up to a numerical, Hamiltonian- independent, coefficient, see Eq. (j6.4ip . The superalgebraic 
structure also detects the difference of the non-coinciding translation parameters. 



6.6 Generic case of n = 2 exact isospectrality 

Consider a generic case of exact isospectrality characterized by the relations ki = k'^, K2 = ^2, 
Ti 7^ r(, r2 7^ T2, see Fig. [5b. The second order intertwining operator possesses then two 
distinct reductions, ()6.2ip and ()6.25p . in which it is necessary to put in addition, respectively, 
K,2 = ^2 and Ki = k'i- The existence of the two third order intertwining operators means that a 
generic isospectral case is described by the basic intertwining operators of the orders 2 and 3, to 
which the intertwining operators and I4 are reduced. To see this, we note that X^ and X^ are 
the third order operators with the same coefficient —1 before the leading derivative term. Then 
the difference of these two operators has to be an intertwining differential operator of the second 
order. This implies that the coefficient before the leading second order derivative term in the latter 
should be a constant. Taking into account that A2 = —A\ + w and B2 = —B\ + w, and employing 
relations (|4.8p and (|6.16p . we find 

Xi - Xi = (Ci - C2)G2 + {kI - kI)Xi , (6.43) 

where Ci = KicothKi(ri — t{), and C2 = K2CothK2(T2 — T2), 

G2 = - ^ ~ ^^'Zx'^ ~Ix ~^ ~'~ ^'^^ K2 (a; + r2) + w'k2 coth K2{x + T2) + , (6.44) 

Xi = ^ + iw-w')+Ci. (6.45) 

In (j6.44p and (j6.45p w corresponds to the function (j6.3p . and w' is the same function but with Tj 
changed for rj, j = 1,2. From (j6.43p it follows immediately that the case Ci = C2 is special, and 
we shall consider it in the next subsection. So, till the end of this subsection we suppose that 

Ci / C2 . (6.46) 

We obtain then the second order intertwining operator 

Y2 = ^f-^ = G2+ "-^Xi . (6.47) 

Li — L2 Ci — C2 



The operator (|6.47p intertwines H'2 and H2, ^2-f^2 ~ H2Y2, and satisfies the relation = ^2'' 

'2 



where 1^' corresponds to Y2 with the interchanged translation parameters Tj and r', j = 1,2. Y2 
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generates the intertwining relation in the reverse direction. Operator Y2, and any of two third 
order operators, or , play now a role of independent intertwining generators. It is more 
convenient, however, to take a linear combination 

Xs = ^-^y^f^ , (6.48) 

different from that in (j6.43p . as a third order intertwining generator to be independent from Y2. 
Using Eqs. (|6.21|) and (|6.25|) . we find that the generic intertwining operators X^ and I4 are reduced 
here as follows, 

(Ci-C2)X5 = {{Ci-C2)H2+Ci4-C2kI)Xs + {4-^4)CiC2Y2, (6.49) 

{Ci-C2)Yi = {t4-Kl)Xs+{{Ci-C2)H2+CiKl-C2Kl)Y2. (6.50) 

Proceeding from the relations (|6.49p . ()6.50p and the relations, presented in the Appendix, which 
correspond to the products of operators X^, I4 and Z5 with the imposed isospectrality relations 
Kj = Kj, j = 1, 2, one can find all the products of the irreducible intertwining operators Y2, X-^, Y2, 

X^, and Lax operators Z5 and Z'^. With these, one can compute the superalgebra generated by the 
second order, Q2-a, and third order, S2-a, supercharges constructed in terms of Y2 and X3 following 
the same rules as we used before, and by the fifth order bosonic integrals 'P2-a- There is another, 
more simple way to compute the superalgebra. Having in mind that fermionic supercharges are 
matrix differentials operators of orders 2 and 3, the alternative form of superalgebra is generated 
by taking a linear combination of them, = CiQ,2-a + S2-a and = C2Q2;a + S2-^a-, constructed 
from X^ and X^ in correspondence with relations (j6.43p and (j6.48p . 

^1=1 j^ABt ) ' ^2 = ^^3^^ ■ (6-51) 

Modifying further the notations, Fa^^ = F^, Fji'^ = F^ , and using the product relations of the 
operators X^'^ , their conjugate, X^'^\ and Lax operators and Zg, see Appendix, we present 
nonzero superalgebraic relations in a compact form, 

{Fj>\F^y = 2Sabhjhihj + 2eabe'^ACVi, (6.52) 

[P2, Fi^)] = 1^ {{-ly h^h2h^2F!,^^ + ^^%jhM'^) . (6.53) 

Here hi = I-L2 + ^ij = ^2 + CiCj, hj = 1, 2, AC = C2 — Ci, and no summation in the indexes i 
and j is implied in the right hand sides. 

Again, the integral Vi = V2-i transmutes here into the bosonic central charge, and the structure 
coefficients depend on both relative translation parameters via Ci and €2- 

The nonzero superalgebraic relations for the third, S2-ai and second, Q2-a-, order supercharges 
and bosonic integrals ■P2;2 can now easily be obtained from (j6.52p and (16.53P by employing the 

relations Q2-a = (-^a^^ — F^i'^)/ AC, S2-a = (C2Fi^^ — CiFji'^) / AC. The superalgebra has the same 
structure (j4.9p . (j4.10p . (j4.1ip as for the n = 1 isospectral case, but with Hamiltonian-dependent 
coefficients of a more complicated form. 
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6.7 Special case of isospectrality with Ci = C2 

Let us consider the special case of isospectrality characterized by the relation 

Ci = C2 . (6.54) 

Equation (j6.54p means that there is a special correlation between relative displacements ri — r{ and 
T2 — T2 and scaling parameters, ki coth ki(t — t[) = K2 coth K2{t2 — T2). In correspondence with this 
relation, we may take a two gap system H2 defined by arbitrary parameters > ki, and arbitrary, 
but finite, ri and T2. Particularly, we can choose the two-gap Poschl- Teller system defined by the 
relations K2 = 2ki and ri = T2- The partner Hamiltonian H2 is given then by the same scaling 
parameters, the finite parameter T2 may be chosen in an arbitrary way with the only restriction 
r2 / T2, while t( is fixed uniquely, r{ = ri — -^arccoth (^ cothK2(T2 — Tg)). 

As a consequence of relation (I6.43P , here a difference — X^ reduces to the first order inter- 
twining operator ()6.45p . which satisfies a relation x\{k,t,t') = —Xi{k,t',t) = —X[. Moreover, 
we will show below that each of the third order intertwining operators X^ and X^ is reducible, 
and so, here the irreducible intertwining operators are Xi and I4. 

As the intertwining generator Xi is the first order differential operator, let us define a superpo- 
tential W by means of 

X^ = ^ + W, W = w-w' + Ci. (6.55) 
ax 

In accordance with relations ([Sj]), KWh . we have W"^ + W' = V2+ Cf, W'^-W = + Cf, 

and then 

XiXl = H2 + Cf, X^Xi = H'2 + Cl , (6.56) 

and X1H2 = H2X1, xIH2 = H2XI. The first order intertwining operator Xi has a form similar 
to that of the operator Xi in one-gap isospectral case. The superpotential W{x) plays here a role 
of the gap function A mentioned there in the context of its relation to the Bogoliubov-de Gennes 
system. 

Operator Xi together with the first order operators X^ = Xi{ki,ti,t[), X^ = Xi{Hi2,T2,T2) 
satisfies in addition the identities 

XfA'^ = aIXi , A2Xf = XiA'2 , X^B'^ = bIXi , B2X{' = XiB'2 . (6.57) 

Let us stress that like ()6.56p . these relations are valid only in the special isospectral case ()6.54p . 
Employing them, we find that the third order intertwining generators X^ and X^ are reducible, 

xi = {H2 + , X^ = {H2 + kI)Xi . (6.58) 

As a consequence, the fifth order generic intertwining operator also is reducible, X^ = {H2 + 
kI){H2 + 4)Xi-CiY^. 

Applying the product relations ()A.32p - ()A.35P collected in Appendix, we can compute the su- 
peralgebra generated by the fermionic supercharges S2;a constructed from Xi and x|, by the 
supercharges Q2-a composed from and vj, and by the bosonic integrals 'P2-a constructed from 



Lax operators Z5 and Z'^^. The nontrivial (anti) commutations relations are 

{SaJb} = 26abhc, , {Qa, Qb} = 25abh\hl , (6.59) 

{Sa, Qb} = 25abCihih2 + 2eabVi , (6.60) 

[V2,Sa] = 2i{hc,Qa - Cihih2Sa) , [^2, Qa] = 2i/li /l2 (Ci " /ll/l25a) , (6.61) 



where hi = 7^2 + j i = 1; 2, hc^ = 'H2 + C^, and we omitted the index n = 2 in the integrals. 
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Supercharges S2-ai cl = 1,2, generate a Lie sub-superalgebra of = 2 supersymmetry. Since 
Cf = C2 > 1^2, it corresponds to the spontaneously broken phase. However, a pecuharity of the 
extended system 7^2 is that it has a structure of centrally extended = 4 nonlinear superasymmetry 
with the two additional fourth order supercharges Q2;a; and two bosonic integrals 'P2-a- Again, the 
integral 1^2-1 plays here the role of the central charge. As in a generic isospectral case from the 
previous subsection, the sum of differential orders of the basic irreducible intertwining operators 
equals 5 and coincides with the order of Lax operators. Again, the superalgebra ()6.59p . ()6.60p . 
(j6.6ip has a remarkable similarity with that for the n = 1 isospectral case. 

We conclude that with a chosen subsystem H2, Eq. (I6.54p defines a one-parametric family, in 
which T2, , T2 ^ r2, is a free parameter of the exactly isospectral system H2. Such a family of the 
Schrodinger pairs is described by the supersymmetry with the two first order supercharges, two 
supercharges of order four, and two bosonic integrals of differential order five, one of which is a 
central charge. This generalizes the one-gap self-isospectral case discussed in Section 4 for the case 
of two-gap isospectral, but not self-isospectral, pairs. 



7 Partially broken and exact isospectralities in n-gap systems 

The analysis of partially broken and exact isospectralities can be generalized for n-gap extended 
systems with n > 2. The case n = 2 considered in the previous Section shows that the concrete 
structure of supersymmetry, namely its irreducible generators and coefficients in the superalgebra, 
depends not only on how many scaling parameters coincide, but also on whether they correspond 
to the same or different ordinal numbers of discrete energy levels of subsystems. It also depends 
on relative translation parameters associated with the corresponding coinciding discrete energy 
levels, and may change in the cases when such relative translation parameters turn into zero, or are 
correlated via equalities of the form ()6.54p . Correspondingly, a concrete form of supersymmetric 
structure is rather variable, but the general picture can be summarized as follows. The n-gap 
pair is characterized by two irreducible basic intertwining operators, one of which is a differential 
operator of odd order, while another is of even order. Each n-gap subsystem also is characterized 
by a nontrivial integral to be a differential Lax operator of order 2n + 1. The orders of irreducible 
intertwining operators satisfy the following rules. As we saw, the case of complete isospectrality 
breaking, when all the scaling parameters of one subsystem are different from those of the second 
subsystem, the supersymmetric pair is characterized by intertwining operators, A'2n+i and Y2n, 
of differential orders |X2n+i| = 2n + 1 and |l2n| = 2n. The sum of their differential orders, 
4n + 1, coincides with the order of the composite differential operator of the form When 
any pair of the scaling parameters of the subsystems coincides, the total order of the two basic 
irreducible intertwining operators decreases in such a way that l-^l'^l = \ {Hn)"'~^P\ = 4n — 1. Any 
new coincidence of some new pair of scaling parameters decreases the total order of XY^ in two. 
Finally, in the case of exact isospectrality, when all the n pairs of the scaling parameters coincide, 
we have \XY^ = \Zn\ = (4n + 1) - 2n = 2n + 1. 

As an example, consider a generic case of exact isospectrality for the pair of the refiectionless 
soliton systems, each having three bound states. In this case, the composite operator A3 has 
6 different factorizations in dependence on the order of the free particle non-physical states ipj, 
j = 1, 2, 3, which are used to generate a 3-gap system. For instance, factorization A3 = A^^ a'^^ A^^^ 
corresponds to that described in Section [5l while A3 = A^^^ A^^ A^^ corresponds to alternative 
factorization like that described in Section 16.21 with A\ constructed in terms of the state ip2, 
A^"^ constructed recursively in terms of A^'^ and Tpi, and finally, A^^^ is constructed recursively 
by employing A^\ A^2^ and tp^. In other words, the upper index indicates here the index of a 
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state ipj we use to construct the first order Darboux operator of tlie generation marked by the 

lower index. The factorizations different from the standard one A3 = Af^A'-^^A^^^ correspond 
to permutations of columns in the Wronskian (j2.2p . and in accordance with Eq. (j2.ip . do not 
produce any effect on the final form of the three-gap potential V3. Employing the information 
on intertwining operators of the n = 2 case, we construct three intertwining operators of order 

where ^3^^^ = 

A2X[^^ A'2'^^'^ and x[^^ = A^^^ — A'j^^'^ — Aq^ is the first order operator constructed in accordance 
with Eq. ()4.2p . and Cr = Kr cothKr(Tr — r^), r = 1, 2, 3. The generic intertwining operator of order 
7 reduces as 

Xj = {Hs + ti^.)xP - CrYe , r = 1,2,3. (7.1) 

Taking (^5'^'* — X^^) and (^5^^ — ^5^^); we get two intertwining operators of order 4, and 
Y^^^\ in which the coefficients before leading derivative term /dx^ will be constants. Presenting 
y(i2) and y(23) 

in a normal form, with leading coefficients to be equal to 1, and taking a difference 
of the resulting fourth order differential operators, we get irreducible intertwining operator of order 
3. Taking any one of the obtained two fourth order operators, we identify finally a pair of the 
basic irreducible intertwining operators X3 and of orders 3 and 4. Three identities in ()7.ip 
allow us then, on the one hand, to express the generic intertwining operators Xj and Yq, which 
are reducible here, in terms of X^ and I4 multiplied by certain polynomials in ^^3. On the other 
hand, the same identities ()7.ip indicate that the cases with Ci = C2 and/or C2 = C3 are peculiar. 
Coherently with the analysis of the previous Section, one can expect that in the special case 
Ci = C2 = C3 the basic irreducible intertwining operators are of orders 1 and 6. The analysis of 
this special case requires a separate consideration and we do not present it here, but only note 
that a corresponding isospectral pair is constructed similarly to the case of the n = 2. Namely, the 
scaling, /C3 > /C2 > ^i, and translation, ti, T2 and r3, parameters of the subsystem i?3 are taken 
arbitrarily, the scaling parameters of the partner system //g are the same, and parameter Tg can 
take any finite value restricted by the condition T3 7^ T3 . The relation C2 = C3 defines T2 uniquely in 
terms of the already chosen parameters, and then the equality Ci = C2 fixes uniquely the remaining 
displacement parameter r(. 



8 Discussion and outlook 

A generic super symmetric quantum mechanical system with a 2 x 2 matrix Hamiltonian, whose 
components are intertwined either by first order Darboux or higher order Crum-Darboux differential 
operators, is described by two fermionic supercharges constructed from the intertwining generators. 
The supercharges together with the matrix Hamiltonian generate, respectively, either linear or 
nonlinear N = 2 superalgebra. We studied a special class of reflectionless systems with super- 
partners having the same number n of discrete energy levels in their spectra. Each of super-partner 
potentials describes an n-soliton solution of a nonlinear KdV equation that depends on n scaling 
and n translation parameters, and satisfies corresponding higher stationary equation of the KdV 
hierarchy. Because of the peculiar, finite- gap nature of the composite matrix Hamiltonians, their 
supersymmetric structure, on the one hand, turns out to be more rich in comparison with a generic 
case, and, on the other hand, it experiences essential changes depending on relation between the 
two sets of 2n parameters that characterize the partner n-soliton potentials. 

There are four supercharges, two of which are composed from intertwining generators X2k+i 
and -^2fc+i differential operators of the odd order 2k + 1 < 2n + 1, while two other fermionic 

integrals are constructed from intertwining generators Y21 and of the even order 21 < 2n, such 
that in general case the total order, |X2fc+i| + 1^2/1; of the basic irreducible intertwining operators 
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satisfies a relation 2n + 1 < {2k + 1) + 2/ < 4n + 1. The system also possesses two bosonic 
diagonal matrix integrals composed from nontrivial Lax operators of the n-gap subsystems, ^2n+i 
and ^2n+i' which are differential operators of order 2n + 1. When the two sets of the scaling 
parameters are completely different, we have a complete isospectrality breaking, and the irreducible 
intertwining generators are of the orders 2n + 1 and 2n. In this case X2n+i and l2n intertwine 
the partner Hamiltonians Hn and i/^ via a virtual free particle system. The anticommutator 
between the supercharges of differential orders 2n + 1 and 2n generates the diagonal Lax integral 
Vn;i = — idiag {Z2n+i, -^2n+i) multiplied by the order n polynomial of the matrix Hamiltonian. The 
both bosonic integrals, Vn-i and T'n;2 = cr^Vn]!, commute nontrivially with the supercharges. The 
Hamiltonian Hn of the system plays a role of the multiplicative central charge of the nonlinear 
superalgebra, whose structure is insensible of the translation parameters of the potentials. In the 
simplest case of n = 1, when the scaling parameters ki and k'i of the partner potentials coincide, a 
kind of a channel for a direct 'tunneling' between the partners is opened, the third order operator 

is substituted for the operator of the first order, that intertwines Hi and H[ directly, without 
communication via the virtual free particle system, and bosonic integral Vi-i transmutes into the 
central charge of the superalgebra, whose structure starts to depend on the 'tunneling distance' 
Ti — r{ . In the case n > 1 , each time when any two discrete energy levels of the partner subsystems 
coincide, the basic intertwining operators X and Y undergo a reduction, decreasing their total 
differential order in two, and a dependence on a relative translation parameter associated with 
a pair of coinciding scaling parameters appears in the superalgebraic structure. The details of 
restructuring of supersymmetry generators depend on whether the discrete energy levels of the 
partners of the same or different ordinal numbers do coincide. A structure of supersymmetry also 
suffers abrupt changes in the orders of the basic irreducible intertwining operators, leaving invariant 
their total sum, when the coincidence of translation parameters, associated with the coinciding 
scaling parameters, happens. The supersymmetry also experiences a restructuring for another kind 
of correlation, kj coth (tj —Tj) = Kji coth Kji (tj' — Tji), j 7^ j' i between the translation parameters 
associated with the coinciding pairs of discrete energy levels of the different ordinal numbers, j ^ j' . 

Only in the case of the exact isospectrality of the partners, when all their discrete energy levels 
coincide pairwise, and as a consequence, their transmission scattering amplitudes also coincide, 
the bosonic integral Vna transmutes into the central charge of the superalgebra. In this case the 
total order 2n + 1 of the two basic irreducible intertwining operators X and Y coincides with 
the differential order of bosonic integrals. A particular case of such a situation corresponds to a 
self-isospectral pair of Poschl- Teller systems. From the viewpoint of supersymmetric structure we 
investigated, the self-isospectral Poschl- Teller pairs possess, however, no special properties when 
n > 1, though the special subfamily of the extended systems with exact isospectrality that we 
detected corresponds to a generalization of the n = 1 self-isospectral case. For n > 1, those 
special isospectral pairs with the scaling and translation parameters correlated by means of (n — 1) 
equalities ki coth Ki(ri —t[) = kj coth Kj{Tj — Tj), j = 2, . . . ,n, are described by the basic irreducible 
intertwining generators Xi and l2n- For n > 1, the corresponding isospectral partner potentials 
have a form different from each other, and if one of them is chosen to be a reflectionless Poschl- Teller 
potential with n > 1 bound states, an isospectral partner does not belong to the Poschl- Teller 
hierarchy of potentials. More precisely, we identified and investigated in detail supersymmetric 
structure of such a special pair in the case n = 2, while we provided here only general indications 
that the same happens for n > 2. The special family of the completely isospectral pairs of n- 
soliton systems with n > 2 requires a separate consideration and will be presented elsewhere. The 
property \Xi\ = 1 means that any of the two hermitian supercharges composed from the irreducible 
intertwining generators Xi and x| may be identified as a first order, Dirac type, Bogoliubov-de 
Gennes finite-gap Hamiltonian that belongs to the AKNS integrable hierarchy. 
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Analyzing the changes of supersymmetric structure associated with a coincidence of the scahng 
parameters, or, that is the same, of the bound states energies, we referred to the opening of tunnehng 
channels conventionally. This might correspond nevertheless to real tunneling processes in some 
applications of the exotic supersymmetry, particularly, related to instantons. 

We discussed the exotic supersymmetric structure from the standpoint of a usual Shrodinger 
equation that corresponds to a potential problem for a particle with a constant mass. It would be 
interesting to reinterpret the results from a perspective of a quantum problem for a particle with a 
position-dependent mass [U] having in mind possible applications for condensed matter physics. 

As it was noted, by displacing the coordinate x for a pure imaginary constant, x x + i5, 
our analysis can be generalized for the case of PT-symmetric quantum systems [l5]. Such a 
generalization seems to deserve a special attention as it was proved to be useful for a particular 
case of supersymmetric extensions of reflectionless Poschl- Teller and related systems, that helped 
recently to clarify some peculiarities in the PT-symmetric quantum mechanics [46]. Particularly, 
T^T-symmetric generalization might be useful for applications in quantum optics. 

As we mentioned, one- and two-gap Poschl- Teller systems control the stability of the kink 
solutions in the sine-Gordon, y?^, and other exotic (l-l-l)-dimensional field theoretical models. By 
considering the doublets of these fields with equal or different masses [48l|49], one could expect that 
the studied supersymmetric structure may reveal itself somehow at the level of the symmetries of 
the corresponding kink solutions. 

We discussed supersymmetric structure by choosing the diagonal Pauli matrix as a grading 
operator F. Alternative choices for T related to reflection operators are also possible. They provide 
the identification of the nontrivial integrals of motion as fermionic and bosonic generators in a 
way different from that described here. Particularly, the treatment of Vn;a as odd supercharges is 
possible, see [50l |42l [41] . Supersymmetric structures for alternative choices of T can be computed 
by employing the product relations of the intertwining generators and Lax operators collected in 
Appendix. The alternative choices were useful for identification of the hidden supersymmetric 
structure in the systems described by the first order Bogoliubov-de Gennes Hamiltonian, particu- 
larly, in those associated with the Schrodinger n = 1 isospectral pair considered here |4H I22j. In 
this direction, it seems to be interesting to apply the results on a special case of the two-soliton 
pairs with exact isospectrality studied in Section 6.7 to the physics related to the Gross-Neveu 
model. 

Finally, it would be interesting to generalize our analysis for finite-gap periodic systems, which 
also find many interesting applications in physics [2T1 EH [2H Ell E2]. In that case it seems to be 
natural to restrict the consideration by the isospectral pairs. 
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Appendix 

Here we collect the products of the intertwining operators and Lax operators to be necessary for 
computing the concrete super algebraic relations. 

In one-gap non-isospectral case, ki 7^ k^, the basic products of intertwining operators and Lax 



25 



integrals are 

XsXl = Hi{Hi + Kl){Hi + Kf), XIX3 = H[{H[ + 4){H[ + Kf), (A.l) 

Y^Y^ = {H, + kI){H, + ), Y^Y2 = {H[ + kI){H[ + k'^), (A.2) 

X^Y^ = -Y2XI = {Hi + k'^)Z:,, Y^X^ = -XIY2 = {H[ + kI)Z'^, (A.3) 

Z3X3 = -Hi{Hi + kI)Y2, XgZ^ = -Hi{Hi + ^2)y2, (A.4) 

^3^2 = (^1 + «?)X3, y2^3 = (^1 + , (A.5) 

Z34 = -Z| = H,{H, + Z'.ZI = -Z'i = H[{H[ + Kf )2. (A.6) 

The products X^Z^, Z'r^X^, 1^2^ .^3 ^^^'^ -^'^^2 ^'^^ obtained by Hermitian conjugation of (IA.4P and 
()A.5p . They are given by expressions of the same form but multiphed by —1 because of the property 
Z3 = — Z3, and with substitutions Hi ^ H[, X^ ^ x\ and I2 Y2 . Relations (|A.6p are needed 
for computing of the superalgebraic structures in the case of alternative choices of the grading 
operator. 

In the n = 1 isospectral case ki = k'^, because of reduction (|4.ip . some relations are changed 

for 

XiXl = Hi + C\ xlXi = H[ + C2, (A.7) 
XiY^ = Z3 + C{Hi + 4), Y2X\ = -Z3 + C{Hi + 4), (A.8) 
Z3X1 = XiZ!, = C{Hi + kI)Xi - {Hi + C^)Y2 , (A.9) 
^3^2 = ^2^3 = (^1 + ^DiiHi + >il)Xi - CY2) . (A.IO) 

The products xlz^ = Z'^x\ and Y2Z3 = Z'^Y2 are obtained by Hermtian conjugation of (IA.9P and 
(jA.lOp as in the non-isospectral case. 

For a pair of ?i-gap systems with complete isospectrality breaking the basic products are 

Y2nYl = FnK , X2„+iXt„+i = H^F^K , (ATI) 

X2„+iy^„ = -y2n4«+l = lPn^2n+l , (A.12) 

^2n+1^2n = IPn^2n+l ) ^2n^2n+l = ^n^2n+l , (A. 13) 

Z2n+lX2n+l = —Hn^nY2n , -'^2n+1^2n+l = —Hn^'nY2n , (A. 14) 

ZL+l = -Hn^n, (A. 15) 

where P„ = U?=iiHn + «?), K = nr=i(-^n + '^f )• Other products of the type xl^^^Y2n etc. 
are obtained from these ones via the change kj -h- k'j, Tj o rj with taking into account that 

■^In+l = ~^2n+l; ^2n ^ ^2 n ^""^ ^In+l = ~^ 2n+l- 

For three cases ()6.19p . ()6.22p . and ()6.22p of n = 2 pairs with partial isospectrality breaking, 
the basic product relations are obtained from IjA.lip -l jA.lSp by taking into account the reduction 
relations (|6.2ip . (|6.25p and (j6.28p . The latter are presented in the unified form X^ = hf^-X^^ —CiY^i, 
and then for each of three cases, distinguished by the index / = 1,2,3 for (I6.19p . (j6.22p . (j6.26p . 
respectively, we have 

xix'^ = hc,K^^K,^, Y,Yl = hlK^K,^, (A.16) 
XlYl = h,,^{Z, + CiK^K^), Y^s^ = V^(-Z5 + CiK^K^), (A.17) 
Z5Yi = K^K^iK^xi - CiYi), Y^Z'^ = K^K,^{K^xi - CiY^), (A. 18) 

Z^Xi = K,{CiKM - hc,Y^)^ 4^5 = K'XCiKM - hc,YA), (A.19) 
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where ha = H2 + a^, a = Hi, Hd, k'^, Q. 

The n = 2 partial isospectrahty breaking case ki = K2 7^ ti = t[, shown on Fig. H^i, is 



characterized by the fohowing basic products of the intertwining and Lax operators, 

F/y/^ = K,K,^ , x,xl = H2hlK,K> , (A.20) 

Xsf/^ = -1^2^X5^ = h,, Z5 , y/^Xs = -x]y/ = /i,. z', , (A.21) 

Z5X5 = -H2hl^K,Y^'' , X^Z'^ = -H2hl^K.Y^^ , (A.22) 

Zsy/ = h^.X^ , %^Z'^ = K,X^ . (A.23) 

For the n = 2 isospectral case with a common n = 1 virtual system, when ki = k'^, K2 = ^^2' 
Ti = r{, T2 / Tg, the basic products are 

Y2^%^' = hl,, Y2^%^ = h'l, XiX^^ = hcX^, X^'xi = h'c,h'l, (A.24) 

XiY^^ = Z^ + C2K, , ^2^X3^^ = -Z5 + C2/iK, , (A.25) 

X^%^ = -Z', + C2K^h'^^, Yi^Xi = Z',+C2KX.^ (A.26) 



Z5y2^ = y2^z^ = hl^xi - C2K,K,Y2^, Y^'^z^ = z',y^^ = C2K,K^y^^ - h%x^\ (A.27) 
z,xi = #4 = C2K.K,Xi - hc,hlj2^, x^^z, = z'^x^^ = h'c,h'l%^^ - C2KX.x^^- 

(A.28) 

Here /i^j = i^2+S^, h'^^ = H'^+kI i = 1,2, hc^ = H2+C2, h'^^ = H2+C2, and C2 = ^2 coth K2(T2-r2). 
Relations for the same isospectral case but with T2 = Tg, ti 7^ r( are obtained from these ones by 
interchanging j4 f-)- i?, ki o K2, ti o r2, r( -H' T2 and by, correspondingly, changing C2 ^ Ci. 

In generic n = 2 isospectral case, ki = k'^, K2 = 1^21 ti 7^ '''2 7^ T2, denoting X^^ = X^ and 
= Xf , we have 

X^^X^'^^ = hihjhj - id - Cj)Z5 , = Kh'.h',. + {C^ - Cj)Z'^ , (A.29) 

Z^X^'^ = Xfz'^ = [{-iyhih2hi2Xf + e'^h,,h]xi''^) , (A.30) 

X«tz5 = Z^X»t = ^ [i-iyh[h',h[,X^'^^ + 6^^M./.f X(^)t) (A.31) 

where hi = H2 + , h'^ = H2 + nf, hij = H2 + CiCj, h'-j = H2 + CiCj, AC = C2 — Ci, and no 
summation in i,j = 1,2 is implied on the right hand sides. 
For the n = 2 special isospectral case Ci = C2, 

XiXl = hc,, Xlx^h'c^, Y^Yl = hlX^, YlY^ = h'lh'l, (A.32) 

XiYI = Z5 + Cih^.h^, , y4X| = -Z5 + Cih^.h^, , (A.33) 

Z5X1 = XiZ'^ = Cih^,h,,Xi - hc.YA , X\Z^ = Z'^X\ = h'cjl - Cih'^X^xl , (A.34) 

Y^Z'^ = Z5Y4 = h^,h^,{h^,h^,Xi - CiY^) , Z'^yI = yIz^ = h'^X^iCiYl - h'^^^xl) , (A.35) 

where = F2 + nf, h'^^ = H'^ + k^^ i = 1, 2, h^ =H2 + Cj, h'^^ = + Cj. 
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